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Abstract. The statistics of the resonance widths and the behavior of the survival probability is studied 
in a particular model of quantum chaotic scattering, (a particle in a periodic potential subject to static 
and time-periodic forces) introduced earlier in Ref. [plj6|. The coarse-grained distribution of the resonance 
widths is shown to be in good agreement with the prediction of Random Matrix Theory (RMT). The 
behavior of the survival probability shows, however, some deviation from RMT. 
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1. The abstract Random Matrix Theory (RMT) is 
known as a powerful tool for analyzing complex quantum 
systems. During the last two decades the predictions of the 
hermitian RMT (which is supposed to describe the prop- 
erties of a closed system) were checked for a large number 
of physical models and an understanding of the conditions 
of applicability of RMT was reached. The situation is dif- 
ferent, however, for nonhermitian RMT, which is aimed to 
describe the spectral properties of open systems. Here we 
have quite a few physical models which can serve a suitable 
test for a nonhermitian RMT. Up to our knowledge the 
following systems are mainly under discussion: chaotic 2D 
billiards with attached leads [Q : the kicked rotor with an 
absorbing boundary condition |E|; simplified models of a 
dissociating molecule || ; scattering on graphs [f| ; and the 
Bloch particle affected by static and time periodic forces 
(ID model of a crystal electron in dc and ac fields) 

The latter model has a number of nice features, which 
distinguishes it among the other systems. First, it is sim- 
ple for numerical analysis. Second, it always realizes the 
so-called case of perfect coupling (where the predictions of 
the hermitian and nonhermitian RMT are most different). 
Third, as a physical model it can be and has been studied 
under laboratory conditions |8|. This present paper con- 
tinues our study of the Bloch particle in ac and dc fields 
in relation to RMT. In particular we discuss here the de- 
cay of the probability in the system (probability leakage), 
which is the simplest quantity measured in the laboratory 
experiment. 

2. We briefly recall some of the results of our preceding 
papers After an appropriate rescaling the Hamil- 
tonian of the system of interest can be presented in the 



dimensionless form 



H = 



+ cos x + Fx + F w x cos(wi) 



(1) 



The parameters of the system (Q) are the amplitude of 
the static force F, the amplitude F u and the frequency to 
of the time-periodic force, and the scaled Planck constant 
h which enters the momentum operator p = ~ihd/dx. 
Another, unitary equivalent, form of the Hamiltonian (0) 
reads as 



H = — + cos[x — e cosyuit)] + Fx 



(2) 



When some condition (based on Chirikov's overlap cri- 
terion) on e is satisfied, the classical dynamics of the sys- 
tem (g) is an example of chaotic scattering. In fact, one of 
the main characteristics of the classical chaotic scattering 
is the delay or dwell time 



t = lim 



[r(po -> -Po) - 2P°/F} 



(3) 



(Here t(j>q — > — po) is the time taken by the particle to 
change its initial momentum pq to the opposite one.) Fig- 
ure |l| shows the delay time (||) as a function of the initial 
coordinate x. The fractal behavior is typical for chaotic 
scattering 
/Veprc 

Ref. I^ M t hat the dynamics and spectral properties of the 
system (0) depend crucially on the condition of commen- 
surability between the so-called Bloch period Tb — h/F 
and the period = 2tt/uj of the driving force, 



We proceed with the quantum case. It was shown in 



Tb 



(4) 
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Fig. 1. Fractal structure of the classical delay time (g) as a 
function of the initial coordinate. The system parameters are 
F = 0.3, e = 1.5, and u = 10/6. 



Providing the condition (0) is satisfied, the complex quasi- 
energy spectrum of the system (the resonances) is given 
by the eigenvalues of the following nonunitary matrix 



U, 



( 



MxN OmxM 



(5) 



In Eq. (|5|) Wnxn is the unitary matrix with the coeffi- 
cients 

W n ', n — (n'\ exp(— ikx)W exp(ikx)\n) , (6) 
Tr {p-Ft + hkf 



W = exp < — — 

h JO 



+ V 



dt 



V = cos[x — ecos(wi)] , 

and OmxN, OmxM, Onxai are blocks of zeros. The nonuni- 
tary matrix (||) can be thought of as the truncated [to the 
size (N+M) x (N + M)) unitary matrix of the system evo- 
lution operator over the common period T = qTg = rT u . 
Then the parameter M, which plays the role of the num- 
ber of channels, is identical with the integer q in condition 
(Q) . The additional parameter N measures the number of 
states supported by the chaotic component of the classi- 
cal phase space. (Unimportant for our present aim is the 




Fig. 2. Distribution of the scaled resonance width _T S of the 
nonunitary random matrix U ran - Parameters are N = 41 and 
M = 1. The statistical ensemble involves 5000 matrices. The 
smooth curves are the distribution (^) for M = 1, 2, 3, respec- 
tively. 



quasimomcntum k which can take any value in the interval 
— n/q < k < n/q.) 

The main conjecture made in Ref. j7| is that the spec- 
tral statistics of the system [i.e., the eigenvalues statis- 
tics of the matrix U sys ] is the same as the statistics of the 
eigenvalues of a random matrix U ran of the structure (||) 
but with the matrix Wnxn substituted by a member of 
the Circular Unitary Ensemble (CUE) 



W N 



xN 



A 



NxN 



A 



NxN 



G CUE 



(7) 



In what follows we examine this conjecture in more detail. 

3. First we discuss the spectral statistics of the random 
matrix U ran . There is strong numerical evidence that the 
statistics of the eigenvalues exp(— i£) — exp(—iE — r/2) 
of the matrix U r an is given by the universal distribution 
derived in Ref. (9). (An analytical proof of this result is 
still an open problem. []) In particular, the distribution of 
the scaled resonance width r s = nF j A {A = 2tt/N is the 
mean level spacing) obeys the equation 



n M (r s ) = 



(-irr t 



M rM-1 



(M-l)! dr™ 



sinh r s 



(8) 



As an example Fig. [2] (adopted from paper 0) shows the 
histogram for the distribution of r s for N = 41 and M = 
1. We note that II M (r s ) w M/T^ for r s > 1 and, thus, 
the notion of mean resonance width is not defined. 

Let us discuss now the decay of the probability P(t), 
which is given by the following equations 



P(t) = \c t \ 2 , C t+ 1 = UranCt , |c | = 1 



(9) 



It is obvious that the dynamics of P{t) is determined by 
the spectrum of the system. Thus the behavior of P(t) 

1 After the paper was submitted we learned that this result 
has been proved in Ref. [p.3|. 
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Fig. 3. Decay of the probability for the random matrix model 
(^) for TV = 121 and M = 1, 2, 3. (Statistical ensemble involves 
100 matrices). The dotted and dashed lines correspond to the 
asymptotic (lid). 



suggests an additional test of the eigenvalue statistics |To| . 
The main advantage of studying P(t) is that this quan- 
tity is more easily measured in the laboratory experiments 

s 

The problem of probability decay was considered (al- 
though within a different RMT model) in paper It 
was proved there that the function P(t) has an exponential 
short-term and an algebraic long-term asymptotic. Adopt- 
ing these results to the present model (||) we obtain 



P(t) 



exp( 



M j 

N 



t > t* 



(10) 



where t* ~ N/2 = tt/A is of order of the Heisenberg time. 
The results of a numerical simulation of the dynamics of 
P(t) depicted in Fig. || well support the analytical expres- 
sion (JToJ). 

It is interesting to note that Eq. ( |l0| ) can be obtained 
by using rather simple arguments. In fact, expanding the 
initial vector Co over the set of eigenvectors of the matrix 
U ran and following Ref. jl2| using the diagonal approxi- 
mation we have 



P{t) 



/ n M (r s )ex P (-rt)dr s 

Jo 



(11) 



For the long-term asymptotic only the narrow resonances 
are of importance. Then, substituting 7Tm(^ s ) by its 
asymptotic expression II m (r s ) ~ ^ S M_1 , f s < 1 [see 
Eq. (§)], we obtain P(t) ~ (2t/N) • This power law 
decay takes place only for a "coherent" evolution of the 
initial state. In contrast, the short-term asymptotic of P(t) 
coincides with an "incoherent" evolution, which would 
take place if one used uncorrelated matrices U ran in Eq. (9) 
at each time step. Then, as follows from the structure of 




Fig. 4. Distribution of the scaled resonance width _T 3 of the 
system (jij) for some values of the static force F satisfying the 
commensurability condition (^) . The other system parameters 
are w = 10/6, e = 1.5, and fi = 0.25. The smooth curves are 
the distribution (^) for M = q. 

the matrix U ran , at each time step the state vector de- 
creases its norm by the factor N/ (N + M) and 



P(t) 



N 



N + M 



exp 



M 



(12) 



4. We proceed with the statistics of the resonances 
for the physical model (||). Numerically we find them as 
the eigenvalues of the matrix (|J), where the parameter 
M is identical with the denominator q in the condition of 
commensurability (^). 

Figure || compares the distribution of the scaled res- 
onance widths r s = ttT/A in the system (g) against the 
prediction of RMT given by Eq. (||). It seen that the global 
features of the distribution n(r s ) fit well to the result of 
RMT. (The peak-like peculiarities of ZZ (r s ) are due to 
the resonances associated with the stability islands of the 
classical phase space. In principle these resonances should 
be removed from the analyzed data.) 

We would like to note that to satisfy the condition 
(||) we adjusted the amplitude of the static force F (the 
other system parameters are kept fixed). By changing F 
we actually change the classical properties of the system 
[in particular, the distribution of the classical delay time 
(B)]. Nevertheless, the quantum distribution n(r s ) re- 
mains practically unchanged (see cases F = 0.0095 and 
F = 0.0663, for example) and is defined exclusively by 
the channel number q. This fact clearly demonstrates the 
applicability of RMT for the system under study. 

We come to the problem of the probability decay. In 
our numerical analysis of the system we calculated the 
dynamics of probability P{t) by two different methods. 
The first method utilizes Eq. (m where the random ma- 
trix U ran is substituted by the matrix (M). The second 
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Fig. 5. Decay of the probability of the system (^) for the values 
of F corresponding to r/q = 1, 3/2, 4/3, and 5/4 in Eq. (|). 
The other parameters are the same as in Fig. 4. 



method is the direct numerical simulation of the wave- 
packet dynamics of the system (Q). The latter method has 
the advantage that it allows to study the incommensurate 
case but it is essentially more time consuming. In the com- 
mensurate case Tb/T u = r/q (with relatively small r and 
q) both methods give the same result. 

Figure || shows the behavior of P(t) on a double loga- 
rithmic scale for r/q — 1, 3/2, 4/3, and 5/4. It is seen that 
the survival probability follows asymptotically a power law 
P(t) ~ t~ a . However, the value of a (a w 1, 4/3, 5/3, 2, 



respectively) differs from that predicted by Eq. (10). This 
means that for very small resonance width (not resolved 
on the scale of Fig. |J) the actual distribution n(r s ) de- 
viates from the distribution (|^). For the moment we have 
no explanation for this deviation from RMT. 

It is also worth to note that the algebraic decay dis- 
cussed is actually a transient phenomenon for physical sys- 
tems. The point is that RMT deals with an infinite ensem- 
ble while in practice it is always finite (and often consists 
of a single representative) . For a finite ensemble most nar- 
row resonance exists and, thus, a very far asymptotic is 
again an exponential decay with the increment given by 
the width of this most narrow resonance. 
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